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ABSTRACT 

Aims. We investigate the molecular gas properties of the deeply obscured luminous infrared galaxy NGC 4418. We address the 
excitation of the complex molecule HC3N to determine whether its unusually luminous emission is related to the nature of the buried 
nuclear source. 

Methods. We use IRAM 30m and JCMT observations of rotational and vibrational lines of HC3N to model the excitation of the 
molecule by means of rotational diagrams. 

Results. We report the first confirmed extragalactic detection of vibrational lines of HC3N. We detect 6 different rotational transitions 
ranging from 7=10-9 to 7=30-29 in the ground vibrational state and obtain a tentative detection of the 7=38-37 line. We also detect 7 
rotational transitions of the vibrationally excited states v 6 and Vj, with angular momenta ranging from 7=10-9 to 28-27. The energies 
of the upper states of the observed transitions range from 20 to 850 K. In the optically thin regime, we find that the rotational transitions 
of the vibrational ground state can be fitted for two temperatures, 30 K and 260 K, while the vibrationally excited levels can be fitted 
for a rotational temperature of 90 K and a vibrational temperature of 500 K. In the inner 300 pc of NGC 4418, we estimate a high 
HC3N abundance, of the order of 10~ 7 . 

Conclusions. The excitation of the HC3N molecule responds strongly to the intense radiation field and the presence of warm, dense 
gas and dust at the center of NGC 4418. The intense HC3N line emission is a result of both high abundances and excitation. The 
properties of the HC3N emitting gas are similar to those found for hot cores in Sgr B2, which implies that the nucleus (< 300 pc) of 
NGC 4418 is reminiscent of a hot core. The potential presence of a compact, hot component (T=500 K) is also discussed. 

Key words, galaxies: evolution — galaxies: individual: NGC 4418 — galaxies: starburst — galaxies: active — radio lines: ISM — 
ISM: molecules 



1. Introduction 

Luminous infrared galaxies (LIRGs) are intriguing challenges 
to modern astronomy. They emit most of their radiation in the 
infrared (IR) region of the spectrum in the form of dust thermal 
continuum, and have typical IR luminosities of Ljr >10 10 L q . 
In many of these objects, the central power source responsible 
for the total energy output is buried deep inside the dusty central 
region and has an origin that remains unclear. For instance, 
Spoon et al. (2007); Aalto et al. (2007) suggest that dusty, 
compact LIRGs may represent the early obscured stages of 
either active galactic nuclei (AGNs) or starbursts and thus play 
a fundamental role in galaxy formation and evolution. 
Owing to the large column of intervening material, observations 
at IR or shorter wavelengths only detect the surface of the 
optically thick nuclear regions, where molecular gas column 
densities can reach values that exceed A^Fb) = 10 24 cirT 2 . Thus 
optical emission becomes completely absorbed by the gas and 
dust, unless the source geometry allows the emission to escape. 
These high masses of gas and dust may also cause an AGN to 
become Compton thick and hard X-rays to become absorbed by 
the intervening material, adding to problems in identifying the 
nuclear source. At radio wavelengths, free-free emission from 
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intervening ionized material may also obscure an AGN (e.g., 
Sanders & Mirabel 1996). This leads to great observational 
ambiguities that are reflected by the classification of LIRGs 
varying considerably with the observed frequency band. To 
produce the high central luminosity observed in many LIRGs 
requires either an AGN or a nuclear starburst (or a combination 
of the two) to heat a large volume of dust. Dust temperatures in 
the inner hundred to a few hundred parsec vary among galaxies 
and can range from 20-30 K to in excess of 170 K (e.g., Evans 
et al. 2003; Downes & Eckart 2007). Dissipation of turbulence, 
which acts on the gas phase (e.g., Guesten et al. 1993) may 
also indirectly contribute to the heating of dust. Because of 
the aforementioned large amount of obscuring material, the 
interplay between the possible energy sources remains difficult 
to ascertain by direct investigation. 

Millimetre observations of molecular lines provide a po- 
tentially valuable tool for determining the effects of the central 
power source on the interstellar medium (ISM) of LIRGs and 
indirectly constraining some of their key properties, such as 
gas temperature, density, and chemistry. Chemical models 
(Meijerink et al. 2007) show how X-ray dominated regions 
(XDRs), generally expected in the case of accretion onto a 
central compact object, and photodissociation regions (PDRs), 
generated by large UV fluxes from young stars, leave different 
imprints in the ISM composition. The peculiar molecular 
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chemistry of dense hot cores around nascent stars was also 
described by Viti & Bayet (2008). The gas temperature structure 
is also expected to be different around an AGN compared to a 
starburst. For example, in an XDR model bulk gas temperatures 
can be as high as 200 K, while in a starburst the temperatures 
should be around 20-50 K (Meijerink et al. 2007). 

Molecular lines surveys of LIRGs, focused mainly on high 
density tracers, such as HCN, HCO + , HNC, and CS, have been 
carried out by several groups (e.g., Krips et al. 2008; Gracia- 
Carpio et al. 2008; Baan et al. 2008; Imanishi et al. 2007; Gao & 
Solomon 2004; Aalto et al. 2002). However, even if these studies 
provide unprecedented insights into the ISM of active galaxies, 
the interpretation of their results remains debated and more sen- 
sitive tracers need to be found. 

In surveys of external galaxies (see Sect. 4.5), bright HC3N 
7=10-9 line emission is found in a subset of IR luminous galax- 
ies. The HC3N line emission is a useful tracer of warm and dense 
regions and is extremely sensitive to a strong IR-field because of 
its multiple, mid-IR bending modes. The molecule is often found 
to be abundant in Galactic hot cores (e.g., de Vicente et al. 2000) 
and can be easily destroyed by intense UV and particle radiation. 
Thus HC3N serves as a tracer of the gas (and dust) properties of 
galaxies with intense IR fields but where the dense gas is not too 
exposed to destructive radiation. 

The edge-on, Sa-type galaxy NGC 4418 has one of the high- 
est luminosities in HC3N (relative to HCN) found for an external 
galaxy so far (Aalto et al. 2007; Monje & Aalto 2008). Its un- 
usual HC3N emission was first reported by Aalto et al. (2007). 
The inner region of NGC 4418 is deeply dust-enshrouded 
(Spoon et al. 2001) with mid-IR intensities that are indicative 
of dust temperatures of 85 K (Evans et al. 2003) inside a ra- 
dius of 50 pc. The IR luminosity-to-molecular gas mass ratio is 
high for a non-ULIRG galaxy indicating that intense, compact 
activity is hidden behind the dust. Interferometric observations 
of HCN 7=1-0 by Imanishi et al. (2004) show that the bulk of 
the dense gas is contained in the inner 2", corresponding to a 
region of 300 pc in diameter. NGC 4418 is a FIR-excess galaxy 
with a logarithmic IR-to-radio continuum ratio (q) of 3 (Roussel 
et al. 2003). This excess may be caused by either a young pre- 
supernova starburst or a buried AGN (Aalto et al. 2007; Roussel 
et al. 2003; Imanishi et al. 2004). The luminous HC 3 N signature 
was interpreted as young starburst activity (Aalto et al. 2007). 
However we need to perform more detailed observations and 
modeling of HC3N emission of NGC 4418 to interpret its ori- 
gin accurately .A tentative detection of a mm vibrational HC3N 
line (Aalto et al. 2007) led us to explore the excitation and abun- 
dance of HC3 N in NGC 44 1 8 and search for a model of the dense 
gas properties in the inner 300 pc of the galaxy. Throughout this 
paper we assume that the observed HC3N line emission emerges 
from a region smaller than or equal to the scale of the interfero- 
metric HCN observations by Imanishi et al. (2004). 
In this paper, we report the first confirmed extragalactic detec- 
tions of vibrationally excited HC 3 N in the LIRG NGC 4418 with 
the IRAM 30m telescope. We report in total 13 different transi- 
tions of HC3N - including vibrationally excited levels - allowing 
us to compile a first model of the excitation and abundance of 
HC3N. We also report a tentative detection of the 345 GHz ro- 
tational 7=38-37 line with the JCMT telescope. In Sect. 2, we 
present the observations and results in terms of line intensities. 
In Sect. 3, we present the results in terms of HC3N rotational 
diagrams and, in Sect. 4, we briefly discuss the interpretation of 
our results and future aims. 



2. Observations 

The observations were carried out in December 2007 and 
July 2008 with the IRAM 30m telescope on Pico Veleta, 
Spain. The coordinates of the observed position, taken 
from the NASA/IPAC Extragalactic Database (NED), are 
12h 26m 54.6s, -00d52m 39s (EQ 2000) . Calibration scans were 
taken every five minutes and pointing checked against nearby 
bright pointing sources every two hours. The nominal pointing 
accuracy was about 2", much smaller than the typical full width 
at half maximum of our beam . For both epochs, observations 
were performed in wobbler switching mode, with a throw of 
120" and 0.5 Hz phase, to optimize the baseline quality. The 
backends consisted of two filterbanks of 1 MHz and 4 MHz 
channel width, connected to the 100 GHz and 230 GHz receivers 
respectively. The weather conditions in the 2007 run were ex- 
cellent with precipitable water vapour of 1 mm or less. Even 
although the summer weather during the 2008 observations led 
to higher system temperatures, the poorer atmospheric transmis- 
sion was counter-balanced by deeper integrations, leading to typ- 
ical rms noise levels of 0.2 mK at 20 km/s resolution. 
In the last observing run, the tuning of the IRAM 30m receivers 
was optimized to include the largest possible number of vibra- 
tionally excited transitions, and the rotational transition of the 
vibrational ground state. This allowed us to compare the new 
data with our previous observations of the vibrational ground 
state, to exclude possible calibration or baseline errors. By com- 
paring the two datasets, we conclude that the calibration errors 
are smaller than or equal to 20%. 

Data analysis was performed with CLASS 1 and X-Spec 2 soft- 
ware packages. A first order baseline was subtracted from all 
spectra, which are shown in Fig. 1. The properties of all the 
HC3N observed transitions are reported in Table 1, together with 
main beam efficiencies, transitions probabilities and energies of 
the upper states. 

3. Results 

3.1. Detected lines 

The rest-frame frequencies of the observed transitions were 
taken from the NIST database of Recommended Rest 
Frequencies for Observed Interstellar Molecular Microwave 
Transitions 3 and cross-checked with laboratory and theoretical 
values from Yamada & Creswell (1986). The notation for the 
vibrational quantum numbers is in the form 7„-7; (v 6 , v 7 ), fol- 
lowing Yamada & Creswell (1986). 

The interaction between the bending angular momentum of 
the vibrationally excited states and the rotational angular mo- 
mentum of the molecule, leads to a /-splitting of the levels. Each 
vibrational state is thus split into two levels, labelled e or / de- 
pending on the wavefunction's parity properties. These parity 
labels are also shown in Table 1 . A useful reference for labelling 
of doubled levels in linear molecules is Brown et al. (1975). 

We detected the HC3N rotational transitions 7=10-9, 16- 
15, 17-16, 25-24, 28-27, and 30-29. For the 7=10-9, 17-16, 
25-24, and 28-27 lines, we also detected rotational transitions 
of the V7=l vibrationally excited levels. For the V6=l lines, we 
generally have upper limits, the only clear detection being in the 
7=25-24 band. The observed HC3N transitions were selected 

1 http://iram.fr/IRAMFR/GILDAS/ 

2 http://www.chalmers.se/rss/oso-en/observations/data-reduction- 
software 

3 http://physics.nist.gov/PhysRefData/Micro/Html/contents.html 
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to minimize blends by major species. We excluded the possibil- 
ity of blending by radio recombination lines as well as by other 
molecules (e.g., methanol). The lack of methanol emission is in- 
teresting and will be discussed in an upcoming paper. 

Thanks to the relatively narrow line widths of about 120 
km/s, line confusion is not a problem in most of the spectra. The 
most crowded situation is found in the 7=17-16 band, where 
multiple Gaussian fits were used to separate the v 7 = le and vg=lf 
lines. Since e and / transitions come from the /-doubling of the 
same bending mode, we can assume that they have the same 
intensity. Fixing the v 7 =le intensity to the more accurately de- 
termined V7=l/ value (see Fig. 1), we were able to separate the 
V6=l/ contribution. We note that the fitted value for v^-lf is 
consistent with the measured upper limit for the V6=le line. 

The molecule appears to be highly excited, with very bright 
emission even for large upper angular momenta and energies. 
After realizing this, we checked spectra of other molecules at 
higher frequencies to search for HC3N contamination finding an 
asymmetry in the CO 3-2 spectrum at 345 GHz, taken by Monje 
and Aalto in 2007 with the James Clark Maxwell Telescope in 
Hawaii (see Fig 1, bottom right). We interpreted this asymmetry 
as the HC3N 7=38-37 line, whose fitted values are also reported 
in Table 1 . This tentative detection was not included in our anal- 
ysis. 



3.2. Line widths 

The observed line widths vary from 80 km/s for the 7=28-27, 
V7 = 1 transition to 182 km/s for 7=30-29. If we exclude these 
two extreme cases, the remaining lines show comparable line 
widths, with a mean value of 125 km/s and a standard deviation 
of 12 km/s that is smaller than the spectral resolution, which 
for our data is not less than 20 km/s. Line widths do not show 
significant trends when compared to the energy of the upper 
level of observed transitions and thus suggest that the emission 
is originating in the same dynamical region inside the galaxy. 
The observed line widths are comparable to those inferred from 
interferometric observations of other dense gas tracers, such as 
HCO + and HCN (Imanishi et al. 2004). Unfortunately, these ob- 
servations can barely resolve the 2" emitting region and do not 
offer much information about the spatial distribution and dynam- 
ics of the dense gas. 

It is remarkable that an edge-on galaxy exhibits these narrow 
emission lines. This might imply that the dense gas is mostly out- 
side the very nuclear region. High resolution observations would 
help us to determine the velocity field of the central gas. 

3.3. Population diagram 

The observed line intensities were combined in a population di- 
agram (see Fig. 2) to attempt a first LTE analysis of the emis- 
sion. An extensive description of the population diagram method 
for deriving interstellar gas properties is given by Goldsmith & 
Langer (1999). From the basic theory of molecular emission, it 
can be shown that the population of a certain energy level is 
given by 
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where W is the integrated brightness temperature. This is related 
to the observed integrated TJ listed in Table 1 by 
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where G s and 9b are the angular sizes of the source and the tele- 
scope beam respectively. In our calculations we assume 9 S =2", 
which is the upper limit to the size of the HCN 1-0 emitting 
region found by Imanishi et al. (2004) with interferometric ob- 
servations. We note that we consider source-averaged quantities, 
since the size of the source is much smaller than the beam and a 
beam-averaged approach would lead to an underestimate of the 
true brightness temperature and column density. Einstein's tran- 
sition probability coefficients A u ( were taken from the Cologne 
Database of Molecular Spectroscopy 4 (CDMS). 
At LTE, we have 
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where /V is the total number density of molecules and Z is the 
partition function. For kT > IiBq, i.e., for warm molecular gas, 
we can approximate Z IcT/IiBq, the rotational constant Bq 
being 4.55 GHz for HC3N. 

Plotting the observed values of equation (1) versus the energy 
of the upper level E u , we can find the excitation temperature of 
the molecule by fitting a line to the natural logarithm of Eq. (3). 
The excitation temperature is given by the slope of the fitted line. 

The population diagram for the observed transitions is shown 
in Fig. 2. The line intensities were fitted using a minimum-^- 2 
method. The free parameters of the fit were column density and 
excitation temperature, which could vary, respectively, within 
the ranges 10 12 -10 16 cirT 2 and 10-600 K. The fit includes opacity 
corrections and results in an optically thin scenario, with typical 
optical depths of 0.04 for v=0 and 0.2 for the v 7 =l transitions. 
A discussion about opacity effects for larger molecular column 
densities can be found in Appendix A. 

The excitation of the molecule is described by four temperature 
components. 



3.3.1. The rotational temperatures 

The v=0 levels cannot be fit by a single line and clearly exhibit 
two temperature components at 29 K (for E M <100 K) and 265 
K (for larger E M ). The v 7 =l transitions have an excitation tem- 
perature of 91 K. These are all rotational temperatures referring 
to the excitation of angular momentum states of the molecule. 



3.3.2. The vibrational temperature 

When we compare the populations of different vibrational lev- 
els with the same 7, we determine the vibrational temperature, 
which describes the excitation of the vibrational modes. This can 
be achieved for the 7=25-24 band, for which we have both V6=l 
and V7=l lines. The resulting vibrational temperature is 519 K. 
The accuracy of the excitation temperatures and column den- 
sities estimated by population diagram fitting were tested by 
Monte Carlo simulations. These results assign a confidence limit 
of about 30% and 20% to the derived temperatures and column 
densities, respectively. 



4. Discussion 

In NGC 4418, the emission of HC3N is the strongest ever de- 
tected in an extragalactic source with a HC3N (10-9)/ HCN (1- 
0) line ratio of about 0.4. The excitation of HC3N is also unusual 



T* A dv, 
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http://www.astro.uni-koeln.de/cdms/ 
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Fig. 1. Spectra of HC3N in NGC 4418. On the X axis, the velocity scale is shown in km/s. The intensity is in Kelvin, not corrected for main 
beam efficiency and beam size. The central tuning frequencies for each observed band are listed in Table 1. In the two cases (top right, bottom 
right) where significant blending is present, we show the Gaussian fits as dashed lines. See text in Sect. 3.1 for discussion about line identification. 



Table 2. Excitation properties of HC 3 N bending modes from Wyrowski 
et al. (1999). Notice the large critical densities, see discussion in the 
text. 



Transition 




Kt 




vibration 




[p-m] 


[s- 1 ] 


[cm 4 ] 




V 6 =l-0 


20 


0.15 


3x10" 


CCN bend 


v 7 = l-0 


45 


6xl0- 4 


4x 10 8 


CCC bend 



4.1. Models of the excitation of HC3N 

If we assume that HC3N is optically thin, multiple temperature 
components must contribute to the emission. This is possible in 
a galaxy characterized by steep temperature gradients in its inner 
regions. We first examine how the radiation field may affect the 
excitation of the HC3N molecule, and then discuss the possible 
origin of the temperature components. 



with a large number of vibrationally excited lines detected, sug- 
gesting that the excitation of the molecule is strongly affected by 
radiation. 



4.1.1. IR pumping of HC 3 N 

The ground vibrational state of HC3N can be excited by IR 
radiation at A =45 fim or /l=20 //m, respectively to a v 7 = l or 
V6=l state (see Fig. 3). The excited state can then decay through 
the R-branch (AJ--1) or P-branch (A/=+l) to populate a v=0 
level, with a resulting selection rule A/=2 (Carroll & Goldsmith 
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Table 1. Observed HC 3 N Lines Properties 



Transition 
J«-J; (v 6 ,v 7 ) 


Rest Frequency 
[GHz] 


Tuning Frequency 
[GHz] 


HPBW 
["] 


7* ; — ~. 

[mK km/s] 


1 b 

[K] 


AV C 
[km/s] 


tb 


A e . 
[s- f ] 


K 
[K] 


Ref.s 


10-9 (0,0) 


90.979 


91.114 


27 


711 ± 20 


1.5 


134 


0.77 


5.8e-05 


24.1 


IRAM 07-08 


(le,0) 


91.128 


// 


27 


<70 


- 


- 


0.77 


5.8e-05 


736.6 


IRAM 07-08 


(0,le) 


91.202 


// 


27 


159 ± 17 


0.4 


108 


0.77 


5.8e-05 


342.6 


IRAM 07-08 


(0,lf) 


91.333 


// 


27 


155 + 14 


0.4 


110 


0.77 


5.8e-05 


342.6 


IRAM 07-08 


16-15 (0,0) 


145.561 


145.561 


17 


1700 ± 80 


1.7 


130 


0.70 


2.4e-04 


59.4 


Aalto 2007 


17-16 (0,0) 


154.657 


154.948 


16 


469 ± 52 


0.5 


117 


0.67 


2.9e-04 


67.3 


IRAM 07-08 


(le,0) 


155.032 


// 


16 


<150 


- 


- 


0.67 


2.9e-04 


779.2 


IRAM 07-08 


(0,le) 


155.037 


// 


16 


316 + 59 


0.3 


141 


0.67 


2.9e-04 


385.2 


IRAM 07-08 


(0,11) 


155.259 


// 


16 


316 + 59 


0.3 


141 


0.67 


2.9e-04 


385.3 


IRAM 07-08 


(le,0) 


154.911 


" 


16 


<150 


- 


- 


0.67 


2.9e-04 


779.2 


IRAM 07-08 


25-24 (0,0) 


227.418 


227.757 


11 


2020 ± 140 


1.1 


140 


0.54 


9.3e-04 


143.0 


IRAM 12-07 


(le,0) 


227.793 




11 


480 ± 110 


0.3 


125 


0.54 


9.3e-04 


853.9 


IRAM 12-07 


(0,le) 


227.977 


// 


11 


928 ± 120 


0.6 


125 


0.54 


9.3e-04 


459.9 


IRAM 12-07 


28-27 (0,0) 


254.699 


254.699 


10 


1720 ± 100 


0.9 


141 


0.47 


1.3e-03 


178.6 


IRAM 12-07 


28-27 (0,0) 


254.699 


255.018 


10 


1315 ± 100 


0.7 


159 


0.47 


1.3e-03 


178.6 


IRAM 07-08 


(le,0) 


255.119 


tt 


9.7 


<200 






0.47 


1.3e-03 


889.1 


IRAM 07-08 


(0,le) 


255.324 


ft 


9.7 


327 ± 82 


0.3 


83 


0.47 


1.3e-03 


495.2 


IRAM 07-08 


30-29 (0,0) 


272.884 


272.884 


9 


1300 ± 200 


0.5 


182 


0.44 


1.6e-03 


204.6 


IRAM 07-08 


38-37 (0,0) 


345.609 


345.609 


14 


1930 ± 250 h 


1.8 


120 


0.63 


3.3e-03 


326.0 


JCMT 2007 



a Half power beam width. 

b Brightness temperature, assuming a 2" source size. 

c Full width at half maximum from Gaussian fit. 

d Main beam efficiency. 

e Einstein's transition coefficient. 

f Energy above ground of the upper level. 

g Telescope and epoch. 

h From Gaussian fit to the CO 3-2 spectrum, see Fig 1. 



1981). The level population can be estimated from Boltzmann's 
equation 

«« B (u I g u _i !L 

— = = — e (4) 

n { B ut I + A u( g t 

where ( and u refer respectively to the ground and excited vibra- 
tional state. Assuming g(B( u = g u B u (, the maximum pumping 
rate is given by 

B( U I = ■£ — = P[ U - (5) 
ekT — l 

This has to be compared to the collisional excitation rate 
Cj-i^j = n(Hi)q.i-\^j for the rotational levels of the ground 
vibrational state. The previous formula gives an upper limit for 
the pumping critical density of n cp = Peulq_j-\^>j- 
For n(H 2 ) < n cp the radiative pumping process is more efficient 
than collisions in exciting the molecule. 

The collisional coefficients q.j-\^j of the first 20 rotational 
levels of HC3N, for different gas temperatures, can be found in 
the Leiden Atomic and Molecular Database. We considered an 
impinging IR flux coming from an optically thick dust source 
at 85 K. We used this temperature for determining both pump- 
ing rates and collisional coefficients. Typical pumping rates are 
of the order of 10~ 5 s _1 for both A-20 pm and /l=45 pm across 
all the range of upper state angular momenta. The collisional 
coefficients at 85 K do not show significant variations for dif- 
ferent transitions and have a value of about 10~ 10 cirT 3 s _1 . The 
resulting pumping critical density is n cp - 10 5 cm 4 . The den- 
sity structure of the HC3N emitting gas in NGC 4418 remains 



unknown and therefore we cannot provide a quantitative esti- 
mate of the impact of pumping on the ground state rotational 
levels. However, for a clumpy medium where dense clouds are 
surrounded by low-density diffuse molecular gas, the pumping 
would act directly on the diffuse medium raising the global in- 
tensity and rotational temperature of the emission. 

4.1 .2. The cool, 30 K, and warm 90 K components 

The 90 K rotational temperature of the v 7 bending ladder may 
originate in dense gas within the inner 0.5". The dust tempera- 
ture of 85 K found by Evans et al. (2003) agrees well with the 
rotational temperature, which may be because either the IR ra- 
diation field dominates the excitation, or the gas and dust are in 
thermal equilibrium and the rotational temperature reflects the 
gas kinetic temperature. 

The low temperature, 30 K, component could be coming 
from a more extended phase, at a greater distance from the warm 
dust in the center. The existence of this cool component is in- 
ferred from the lower transition v=0 lines and can be confirmed 
by observations of low-/ lines with, for example, the Effelsberg 
telescope. 

4.1 .3. The hot, 260 K and 500 K, components 

As we can see in Fig. 2, the estimated rotational temperature of 
the high-/ transitions is higher than the 30 K found for the low- 
/s, by almost one order of magnitude. This change in the slope 
of the rotational diagram cannot be explained by opacity effects 
and must be investigated further. Concave population diagrams 
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Fig. 2. Population diagram for all the observed HC3N lines in 
NGC 4418. The 3-cr error bars include uncertainties in the line inten- 
sity and beam size. The solid lines represent the fit of the level pop- 
ulation for v=0 (circles), V6=l (triangles), and vj = \ (diamonds). The 
dashed line fits the v=0, v 6 =l, and Vj=l vibrational states of the 7=25 
level (marked with a star). For the V6=l lines, we note that the only 
detection is for 7=25-24, the others being only upper limits. The vibra- 
tional ground state was fitted with two temperature components. The 
fit includes opacity corrections, as described in Appendix A. Maximum 
opacities are 0.04 for v=0 and 0.2 for v 7 =l. 



are commonly found in Galactic molecular observations and are 
often attributed to multiple temperature layers. Temperatures of 
300 K are typical of Galactic hot cores around sites of massive 
star formation (see Sect. 4.4), the excitation temperature found 
could thus be reflecting the true gas kinetic temperature. This is 
in agreement the detections of the 14 fim absorption of HCN, for 
which Lahuis et al. (2007) find gas temperatures of about 300 K 
in the core of NGC 4418. 

The vibrational temperature of the 7=25-24 lines is about 500 K. 
The bending modes vg and V7 both have critical densities greater 
than 10 8 ctrT 3 (see Table 2) and are thus probably radiatively ex- 
cited. The derived vibrational temperature then reflects the tem- 
perature of a radiation field. 

4.1.4. Radiative excitation of high-J levels 

The high temperature of the high-/ rotations may also be caused 
by radiative excitation, that becomes more and more efficient as 



V7 = l 



P-branch 



J-l_ 
J-2 _ 



V 7 =0 




Fig. 3. Radiative pumping of HC 3 N rotational levels via vibrational IR 
transitions. 



the angular momentum of the upper state increases. The critical 
density of a / — > / - 1 transition can be written as 



«c = 



Ajj-i 



64n 4 3 2 J + 1 



1 



_ -V n 

qjj-\ 3hc 3 27 + 3 qjj-i 



(6) 



where v is the line frequency and /j. is the electric dipole moment 
of the molecule. For n(Hi) < n c , the excitation temperature is 
determined by the radiation field and not the gas kinetic temper- 
ature. The critical densities of the first 20 rotational transitions 
of HC3N are shown in Fig. 4. We see that n c increases from 10 3 
ctrT 3 for 7 M =1 to more than 10 6 ctrT 3 for 7„=20. For /„ > 15, 
we note that the critical density is higher than 10 6 cm 4 . This 
region corresponds to energies E u greater than 100 K, which is 
the critical value at which the population diagram of Fig. 2 ex- 
hibits a change in slope towards higher excitation temperatures. 
For n(H 2 ) < 10 6 , this may imply that the excitation of high-/ 
levels is affected considerably by continuum radiation. 



4.1 .5. Can the 500 K component exist? 

Infrared observations by Evans et al. (2003) report a dust tem- 
perature of 85 K. This value was inferred by comparing the IR 
fluxes at 60 and 100 yum, and can be combined with the total 
IR flux to estimate the size of the emitting region. In their pa- 
per, Evans and collaborators find an optically thick dust source 
of about 0.5" in size, corresponding to 70 pc at a distance of 30 
Mpc. 

Assuming that our 500 K component also represents the tem- 
perature of the IR continuum, the required source size 6 for the 
emitting region can be estimated by the relation L mIR cc 2 T 4 , 
which infers 0=0.01". This corresponds to a linear diameter of 
1 .45 pc. This is an upper limit to the size of the emitting region, 
where it the origin of all the observed IR flux. Our interpretation 
of the molecule's excitation then implies that an extremely com- 
pact radiation source exists in the core of NGC 4418. This may 
correspond to a deeply embedded AGN, heating the surrounding 
high column of dust to the required high temperatures. A hot 500 
K component would peak at 10 fim, right on top of the silicate 
absorption band, which in NGC 4418 is one of the deepest ever 
measured in a luminous IR galaxy (Spoon et al. 2007). This high 
temperature feature might be hidden in the silicate dip and not 
be evident in the SED. Furthermore, Lahuis et al. (2007) find gas 
temperatures from the 14 //m absorption of HCN to be 300 K in 
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the core of NGC 4418. This would require the continuum source 
to be hotter than 300 K gas observed in absorption. 

4.2. Models of the excitation of HC 3 N - optically thick case 

If we abandon the notion that HC3N is optically thin, we show 
in Appendix A that optical depth effects may allow us to fit a 
single temperature (with however significant scatter) to the v=0 
transitions. This lowers the excitation temperature of the high-J 
levels, but leaves us with extreme column densities of HC3N and 
large optical depths that should result in strong lines in the 13 C 
isotopic variants of HC3N. 

Assuming a Galactic C/ 13 C abundance ratio, these lines should 
have intensities comparable to the emission of the main iso- 
topomer. However, measurements of the 3 mm band obtained 
using the new EMIR receiver at IRAM 30m (Costagliola et al., 
in prep.) show that the emission from H 13 CCCN J=10-9 is at 
least 10 times fainter than that from the main 12 C variant. Even 
when attempting to correct for optical depth effects, we still have 
trouble fitting all the points to a single temperature and the col- 
umn densities required (7V(HC3N)> 10 17 cm" 2 ) imply HC3N 
abundances in excess of those of HCN and HCO + , starting to 
approach those of CO. 

4.3. HC3N abundance 

From the population diagram in Fig 2, the V7 component has a 
column density of 8.7xl0 15 cirT 2 , which can be assumed to be 
a lower limit to the total HC3N column. 

The hydrogen column is estimated from 13 CO J=l-0 and 
J=2-l measurements, taken in 2007 and 2008 at IRAM 30m 
(Costagliola et al., in prep.). A basic LTE interpretation of 
13 CO intensities leads to a source-averaged column density 
N( 13 CO)^5.2xl0 16 cnr 2 . Assuming a C/ 13 C ratio of 50, typi- 
cal of starburst galaxies (e.g., Henkel & Mauersberger 1993; 
Mao et al. 2000) and applying a N(CO)/N(H 2 )=10" 4 conver- 
sion (Blake et al. 1987), we obtain a molecular hydrogen col- 
umn N(H2)- 2.6 x 10 22 cirT 2 . This estimate is consistent with 
the value of 7.7 xlO 22 cirT 2 derived by Lahuis et al. (2007) from 
silicate absorption in Spitzer spectra. The resulting HC3N abun- 
dance is then X[HC 3 N]=N(HC3N)/N(H 2 )=3.4xl0- 7 . 
In our calculations we assume that the 13 CO emission is optically 
thin and effectively traces the hydrogen column. The observed 
12 CO/ 13 CO J=l-0 ratio in this galaxy is about 20, implying only 
moderate 12 CO opacities. It is therefore reasonable to assume 
optically thin emission for 13 CO. 

4.4. HC 3 N in the Galaxy 

In our Galaxy, HC3N is often associated with hot cores, i.e., 
dense and warm regions around newborn massive OB stars. 
Here, large dust columns (Ay ^1000) provide shielding against 
UV radiation from the central object. The first interstellar detec- 
tion of vibrationally excited HC3N was by Clark et al. (1976) 
toward the Orion hot core and an extensive survey of vibra- 
tionally excited lines in other Galactic sources was performed 
by Wyrowski et al. (1999). Bright HC3N emission has also been 
found towards circumstellar envelopes and regions around plan- 
etary nebulae (Pardo et al. 2004). 

Our estimate of the HC3N abundance has the same order of 
magnitude as that found by de Vicente et al. (2000) towards hot 
cores in Sgr B2, one of the most active regions of high mass star 
formation in the Milky Way. An enhancement of HC3N abun- 



dance in hot cores was also observed in Orion A (Rodriguez- 
Franco et al. 1998), where the abundance is higher by about one 
order of magnitude than that of the diffuse warm gas, reaching 
values close to 10~ 8 . In the same complex, HC3N strongly anti- 
correlates with HII regions, where it is likely to be destroyed by 
the strong UV radiation. 

The estimated HC3N abundance for NGC 4418 thus resembles 
the properties of Galactic hot cores, where the highest concen- 
tration of the molecule has been observed. It is quite remarkable 
that this appears to be a global property of the galaxy, especially 
considering that 13 CO emission might be more extended than 
HC3N and thus the estimated abundance might be only a lower 
limit to its true value. 



4.5. HC3N in other galaxies 

One reason for selecting NGC 4418 to perform a deeper study 
of its HC3N component was the unusual luminosity of its HC3N 
10-9 line. Combining observational data with data in the lit- 
erature, Lindberg (2009) classified galaxies as HC-iN -luminous 
when the HC 3 N 10-9 line was at least 15 % of the HCN 1-0 line 
intensity. He found that only a few (5) of 30 galaxies in his sam- 
ple fitted this category. NGC 4418 belongs to this exclusive class 
of HC3N-luminous objects, having an HC3N/HCN line intensity 
ratio of 0.4. In a high resolution study of IC 342, Meier & Turner 
(2005) find that HC3N emission correlates well with the 3 mm 
continuum emission, and avoids regions of high UV radiation. 
For example, they find a spatial anticorrelation between HC3N 
and the PDR-tracer C2H. The conclusion is that HC3N traces 
young star-forming regions. On the scale of the main molecu- 
lar complexes, the HC3N/HCN ratio is 0.28, which is still lower 
than the global value for NGC 4418. Wang et al. (2004) report a 
summary of HC3N abundances in the nearby starbursts NGC 253 
and M 82 as well as the Seyfert-2 NGC 4945, where these range 
from X=10~ 90 to 10~ 79 . For these three galaxies, Lindberg lists 
the HC3N/HCN intensity ratio to be less than 0.1 for NGC 253 
and NGC 4945 and for M 82 it is less than 0.2. We note that 
these are central and not global ratios, these three galaxies being 
extended compared to the single-dish beam. 

4.6. Outlook 

To confirm the presence of a 260 K component, we need care- 
fully calibrated submm data of the high-/ transitions. To more 
tightly constrain the excitation and spatial extent of HC3N, high 
resolution interferometric observations are important. A hot, 500 
K, dust component should be identifiable by high resolution IR 
observations. 

In a previous paper we proposed that HC3N could be used as 
an indicator of AGN since HC3N should be destroyed by both 
ions and UV radiation. In a follow-up paper, we will present 
models of the impact of XDRs and PDRs on the HC3N abun- 
dances and discuss the survival and formation of HC3N in these 
environments. We do not discuss in this paper the possibility 
that the hot component is not excited by a radiation source, but 
that the temperature instead reflects the formation and destruc- 
tion processes of HC3N, thus indicating that they must be taken 
into account when modeling the excitation of the molecule. We 
are aware of the fact that a non-LTE analysis would be required 
to properly address the problem of the molecule's excitation. A 
complete analysis of the excitation and the chemistry of HC3N 
will be discussed in a follow-up paper. 
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5. Conclusions 

We have confirmed the first extragalactic detection of vibrational 
lines of HC3N. We have detected 6 different rotational transi- 
tions ranging from 7=10-9 to 7=30-29 in the ground vibrational 
state, plus a tentative detection of the 7=38-37 line. We have 
also detected 7 rotational transitions of the vibrationally excited 
states V6 and vq, with angular momenta ranging from 7=10-9 to 
28-27. In the optically thin regime, we find that the v = tran- 
sitions can be reproduced by models of two temperatures, 29 K 
and 265 K, while the vj lines can be fitted by a temperature of 
91 K. The vibrational temperature, fitted to the J=25-24 transi- 
tions is, 519 K. By allowing the column density to vary between 
the different temperature components, we inferred HC3N col- 
umn densities of 8xl0 14 cm" 2 for the vibrational ground level, 
and 9xl0 15 cm" 2 for the vibrational vj-1 transitions. 
The excitation of the HC3N molecule responds strongly to the 
intense radiation field and the presence of warm, dense gas and 
dust at the center of NGC 4418. The intense HC3N line emission 
is a result of both high abundances and excitation. The HC3N 
excitation and abundances seem similar to those found for hot 
cores in Sgr B2 in the Galactic Center. This implies that the nu- 
cleus of NGC 4418 has properties in common with Galactic hot 
cores. It cannot be excluded that the hot (500 K) component may 
be associated with a buried AGN. 
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Appendix A: Opacity effects 

The effects of a finite optical thickness of the observed lines 
on the population diagram have been investigated by different 
authors, e.g., Goldsmith & Langer (1999) and Wyrowski et al. 
(1999). 

If we were to define C T = t/(1 - e~ T ) as the correction due to 
optical depth, Eq. 1 would read 



N u = 



%nkv 2 W 
hc 3 A, 



C T = N M "C T . 



(A.l) 



lit 



For t << 1, C T - 1 and N u - N' hin , but for r > 1, we have 
Cj — t and the observed line intensities lead to an underestimate 
of the derived column density if no correction is applied. 
For a linear molecule, the optical depth of a 7 — > 7 - 1 transition 
is given by 



87rV N 
3h ZA^ 



Je- aJ(J+1> (e 2aJ -l), 



(A.2) 



where \i is the dipole moment of the molecule and a = hBo/kT ex , 
where Bq is the rotational constant. If kT ex » IiBq, the partition 
function can be approximated by Z ^ a" 1 (Goldsmith & Langer 
1999). This is often the case for HC3N in interstellar space, since 
Bo=4.55 GHz and the previous condition becomes 7 ex >>0.2 K. 
The lack of detected emission from 13 C-variants of HC3N (see 
Sect. 4.2) prevents us from obtaining a direct estimate of the 
optical depth. We therefore used Eqs. A.l and A.2 to fit the 
observed intensities via a minimum-^ 2 method. We focused in 
particular on the v=0 lines, to find out whether the change in 
slope could be due to optical depth effects rather than a high- 
temperature population. The result of the fit is shown in Fig. 
A.l. The free parameters were excitation temperature, column 
density, and source size, which were varied across a wide range 
of values. 

We obtain two possible solutions, one at N=8.3xl0 14 cirT 2 
and one at N> 10 17 cm 2 . Both panels in Fig. A.l show that a 
single component does not fit the observed intensities. When we 
reach the optically thick regime (around 10 16 -10 17 cm" 2 ) the fit 
improves very slowly with increasing column density. However, 
even for extremely large columns (e.g. 10 25 cm" 2 ), the normal- 
ized^- 2 is still high (around 20). Thus the change in slope of the 
population diagram cannot be explained by optical depth effects. 

An estimate of the optical depth of the observed lines can be 
obtained by comparing their brightness temperature (T/, in Table 
1) to the excitation temperature. In general, for a spectral line we 
have Tb - (T ex - Tb)(1 - e~ T ). Neglecting the contribution from 
the background T B , the optical depth is given by t = -ln{\ - 
Tb/Tex)- Considering an excitation temperature of 30 K for the 
first three v=0 lines in Table 1, we derive r < 0.06, which is 
consistent with our optically thin scenario. 
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Fig. A.l. Fit of HC3N v=0 intensities including optical depth effects. 
Top: Fit for HC 3 N column density smaller than 10 16 cnT 2 (normalized 
^ 2 =80). Bottom: Fit for HC3N column density larger than 10 16 cm -2 
(normalized ^ 2 =50). In both fits, the source size 6 could vary between 
0.1 and 2 arcseconds. The dashed line represents the optical depth for 
the different transitions, in arbitrary units. 



